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Abstract 

Combination of two basic types of synchronization, anticipated and isochronous synchronization, 
is investigated numerically in coupled semiconductor lasers. Due to the combination, a synchro- 
nization of good quality can be obtained. We study the dependence of the lag time between two 
lasers and the synchronization quality on the converse coupling retardation time Tc2i ■ When Tc2i is 
close to the difference of external cavity round trip time r and coupling retardation time Tc12, the 
combination of anticipated and isochronous synchronization may produce a better synchronization, 
with a lag time proportional to Tc2i- When Tc21 is largely different from r — Tc12, the combination 
is noneffective and even negative in some cases, with a lag time independent of Tc2i- 
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I. INTRODUCTION 



In last few decades, much attention has been paid to chaotic synchronization because of 
its potential applications in wide variety of fields, especially in communications 
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, |l2| , where the intensity output of the laser exhibits 
irregular dropout events (sisyphus effect, see [l^)- Sudden reduction in output (a dropout) 



is followed by a gradual recovery and t 
interval between dropouts is irregular 
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: len another dropout comes 

isl ligj 0, Izil^^cl two lasers without interaction to 
each other should exhibit severally irregular and uncorrelated dropouts. When there is a 
coupling between two lasers by injecting part output of one laser (laserl) into the other 
(laser2), the two lasers may produce correlated outputs, particularly in dropout events 22 1. 
In other words Dropouts of laserl is imaged into the output of laser2, thus dropout events 
of two lasers appear in the same pace, and synchronization between two lasers is achieved in 
this way. Such synchronization doesn't mean always exact equality of two lasers' outputs, 
but can be considered as a chaos control of laser2's LFF behavior by laserl, in spite that 
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complete synchronization can be obtained in special situations 

Through a simple analysis of the rate equations for the unidirectiona 
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system, two types of synchronization has been identified 
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coupling lasers 



Anticipated Synchronization{AS) : /2(t — r) = Ii(t — r^u) 
Isochronous Synchronization{I S) : l2{t) = h{t ~ '^cu) 



(1) 
(2) 



System may choose one type of synchronization to exhibit, and the other synchronization 
behavior is hidden. Which type is chosen is determined by the competition between two 
types. The winner is shown and the loser is hidden. Transition from one type to the other 
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may occur when operating parameters are changed 
only competition but also cooperation or combination between two types of synchroniza- 
tion, especially when the problem is extended to the field of complex network with delayed 



feedback and coupling j3f 



In this paper, after some brief introduction to the two 



basic types of synchronization, the synchronization combination phenomenon is discussed. 
To our best knowledge, this is the first paper to discuss the combination of the two basic 
types of chaotic synchronization in coupled semiconductor lasers system. 



2 



^aser 2 ^ 



7 

M2 

unidirectional Coupling Lasers System 



FIG. 1: Schematic representation of unidirectional coupling lasers system, where two basic types 
of synchronization are shown typically. 

II. ANTICIPATED AND ISOCHRONOUS SYNCHRONIZATION 

The two basic types of synchronization are reviewed in this section. 

The unidirectional coupling lasers system is shown in Fig.l. Part output of laser 1 is 
injected back as a feedback by a mirror Ml. Another part is injected into laser2 via a beam 
splitter BSl, an optical isolator Oil (ensure there is no light from laser2 to enter laser 1 and 
alter the dynamics of laserl) and another mirror M2. No feedback is used in laser2. Photo 
diodes (PDl, 2) is used to detect outputs of two lasers respectively. Feedback rate is labelled 
by 7, ?7i2 is coupling strength from laserl to laser2. 

Numerical simulation is performed by LK equations for the complex electric fields E and 
normalized carrier densities js^. 

' A;(l + za)[G'i-l]Ei(t)+7iEi(t-r)e(-'"i")+/36(t) (3) 



dt 

dNi i-Ni-Gi\Ei 



dt T„, 
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Where subscripts 1 and 2 denote Laserl and Laser2 respectively. The second term in Eq. 
(1) corresponds to the feedback in Laserl, and the second term in Eq. (3) corresponds to 
the coupling from laserl to laser2. k is the cavity loss, a the linewidth enhancement factor, 
G = N/{1 + e\E\^) is the optical gain, e is the gain saturation coefficient, lu is the optical 
frequency without feedback, Au = uj2 — uji is the frequency detuning, ^ is independent 
complex Gaussian white noise, and (3 measures the noise strength, j is the normalized 
injection current, and r„ is the carrier lifetime. For simplicity, we may choose /?i = /?2 = 0. 
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FIG. 2: Time traces of two unidirectional coupling lasers operating in the LFF regime, typically 
showing anticipated synchronization with parameters: ryi2 = 7 = 8ns~^, r = 7ns, Tc12 = 4ns, 
j = 1.003, a = 5, k = 500ns~^, t„ = Ins. Vertically shifted low-pass-filtered ones are plotted in 
solid lines to emulate the experimental detection and exhibit the dropout events clearly [l^ . 

Two types of synchronization are shown in following two cases respectively. 

In the first typical AS (anticipated synchronization) is shown with parameters: 

V12 = 7 = 8ns~^, T = 7ns (r is external cavity round trip time in laserl), Tc12 = 4ns (rci2 
is the time for light to travel from laserl to laser2). The intensity outputs of two lasers 
are plotted in Fig. 2. The two lasers are both operating in LFF regime, where sudden 
output reduction (a dropout) appears, followed by a gradual recovery, and then another 
dropout comes. Intervals between dropouts is irregular. To remove the fluctuations in high 
frequencies so as to make the dropout events clearer, the low-pass-fiitter diagrams are also 
plotted and shifted upward. The dropout events is labelled by short dotted lines. 

As clearly shown in Fig. 2, the two lasers undergo dropouts in the same pace due to the 
coupling and synchronization is achieved. In addition, 3ns before every dropout in laserl, 
there is always a corresponding dropout in laser2, as if laser2 can predict the future state 
of laserl in spite that laser2's chaotic dynamics is driven by laserl. This is a typical AS, 
satisfying Eq. (1). 

In Fig. 3 correlation function C between two lasers' outputs as a function of shift time 
(ts) is also plotted. C is defined as follows: 

{[hit + r,) - {h)][h{t) ~ {h)]) 



(7) 



At Ts = 3ns there is a peak. This peak indicates that a strong correlation can be obtained 
when laser2 is shifted backward by 3ns. Obviously this peak corresponds to the AS in Fig. 
2. 
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FIG. 3: plot of correlation function C as a function of shift time r^, corresponding to the typical 
anticipated synchronization in Fig. 2. 
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FIG. 4: Time traces of two unidirectional coupling lasers typically showing isochronous synchro- 
nization with parameters: r/12 = 5ns~^, other parameters are the same as in Fig. 2. 

Besides the right peak, there is another peak at ^ —4ns. (It is noted that Tc12 = 4ns 
in our simulation. This peak indicates that shifting laser2 forward by r^u also can produce 
a strong correlation. This left peak satisfies Eq. (1) and corresponds to IS. Although only 
AS is shown in Fig. 2, two-peaks configuration in correlation plot implies the existence of 
IS. Only because the right peak is higher (indicating that the anticipated synchronization 
quality is better), AS is the winner in the competition with IS, and IS behavior is hidden. 

In the second case, a typical IS (isochronous synchronization) is shown with parameters: 
rii2 = Sns"^, 7 = 5ns^^, r = 7ns, = 4ns, namely only 7 is decreased a little. Outputs 
of two lasers and vertically-shifted low-pass-filtered diagrams are plotted in Fig. 4, ranging 
from 200ns to 450ns. Every dropout in laserl leads the corresponding dropout by about 4ns 
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FIG. 5: plot of correlation function C as a function of shift time Tg, corresponding to the typical 
isochronous synchronization in Fig. 4. 

(We note Tc12 = 4ns). So there is a typical IS between the two lasers, satisfying Eq. (2). 
This synchronization is induced by the fact that the unidirectional coupling is sufficient for 
the laserl to drive the dynamics of the laser2 and lead to a locking-state phenomenon. 

Fig. 5 is the correlation plot. It is very easy to identify the two peaks. The left peak 
corresponds to IS and the right one corresponds to AS. Only because the left peak is higher 
than the right one (indicating that the isochronous synchronization quality is better), AS is 
the loser in the competition with IS and has been hidden. Hence it is clearly seen in Fig. 4 
that only IS is shown. 

According to above discussion, which type of synchronization is shown, AS or IS, is 
determined by their competition. The winner's behavior is shown, while the loser still exists 
in the system but has been hidden. 

III. SYNCHRONIZATION COMBINATION 

AS and IS not only compete but also cooperate with each other, i.e. their combination. In 
this section such combination is discussed, the model under our following discussion is almost 
the same as the second case in the previous section except that a converse coupling with 
coupling strength 7721 is used as shown in Fig. 6 with parameters: 1721 = V12 = 8ns~^, 7 = 
bns"^. Two optical isolators are used to make sure that both couplings are unidirectional. 
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FIG. 6: Schematic representation of bidirectional coupling lasers system, where the combination 
of anticipated and isochronous synchronization is investigated. The coupling from laserl to laser2 
travels via BSl, Oil and M2, while the converse coupling from laser2 to laserl travels via BS3, 
012 and BS2. 

then Tc2i 7^ Tci2, which is helpful in following analysis. Tc2i is the time taken for the light to 
travel from laser2 back to laserl. 

Then the right correlation peak appears because of a combination of two factors: 

(1) Before adding the converse coupling, there has been an AS: l2{t — r) = — rci2), 
where laser2 lags laserl by the time = r — Tc12- 

(2) After adding the converse coupling, there is also an IS due to the converse coupling 
from laser2 to laserl: Ii{t) = l2{t — Tc2i), where laser2 lags laserl by the time = Tc2i- 

We may choose: r = 7ns, r^u = 4ns, and Tc2i = Sns, then Tc2i = t — Tcu, and there will 
be a combination of AS with IS at = 3ns. 

The dynamical behavior of bidirectional coupling laserl (with a external feedback) and 
laser2 (solitary) is described by Eq. (8, 9) and Eq. (5, 6) respectively. 



k{l + ia)[Gi - l]Ei{t) + -fiEiit - r)e(-^"^^) + r/2i^2(t - re2i)e[-*('^^"=^^-^"*)] + 



dt 

dNi _ j -Ni- Gi\Ei\ 

dt Tn 



(9) 



where the third term in Eq. (8) is added for the converse coupling. 

Numerical simulation is shown in Figs. 7 and 8. As expected, in Fig. (7) 3ns before every 
dropouts of laserl there has been a dropout in laser2 inevitably. Synchronization between 
two lasers is achieved and Laser2 leads laserl by 3ns. 

That Laser2 can predict the future dynamical behavior of laserl results from the combi- 
nation of two basic types of synchronization behavior, AS and IS. The effect of combination 
is also exhibited by the correlation plot Fig. 8. The right correlation peak at = 3ns, re- 
flecting the effect of combination of AS and IS, is higher. While the left peak at ~ — 4n.s 
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FIG. 7: Time traces of two lasers, showing the effect of the combination of anticipated and 



isochronous synchronization, laser2 is seen to synchronize and lead laser 1. parameters are: 7721 = 



8ns ^, Tc2i = 3ns, and other parameters are the same as in Fig. 4. 
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FIG. 8: plot of correlation function C as a function of shift time Tg, corresponding to the combi- 
nation of anticipated and isochronous synchronization in Fig. 7. 

reflecting IS produced by only the coupling from laserl to laser2 is lower and its correspond- 
ing dynamical behavior is hidden. So only the combination synchronization is shown. 

It seems that the matching condition Tc2i = t — makes the combination of AS and 
IS create a strong correlation and consequent combination synchronization. For further 
discussion about the effect of combination of AS and IS, a more general situation where 
Tc2i 7^ T — Tci2 IS conccmed by moving the branch of the converse coupling so that Tc2i is 
modified. 

In the following discussion, combination synchronization quality Qc and its characteristic 
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FIG. 9: shows the dependence of Tc and Qc on the converse couphng retardation time Tc2i ranged 
from 2.4ns to 3.6ns. Five regions are identified. All the parameters are the same as in Fig. 7 
except for the variable Tc2i- 

time scale Tc are investigated: 

Q, = max(C(r,)), r, > (10) 
C(T,)) = (11) 

where Tc is the lag time between two lasers. 

(1) In AS, laser 1 lags laser2 by the time dependent on the difference of the external cavity 
round trip time of laserl and the coupling retardation time , i.e. t — t^u, independent on the 
converse coupling retardation time Tc2i- When Tc2i is increased, the lag time should almost 
remain unchanged. 

(2) In IS, laserl lags laser2 by the time dependent on the converse coupling retardation 
time Tc2i- When Tc2i is increased, the lag time should be increased proportionally. 

In combination of AS and IS, whether the lag between two lasers (Tc) is changed propor- 
tionally or unchanged when Tc2i is increased? 

Fig. 9 shows the dependence of Tc (a) and combination synchronization quality Qc (b) 
on the converse coupling retardation time rc2i. The interesting results show the range of 
Tc2i from 2.4ns to 3.6ns can be divided into five regions, labelled by I-V. In region I and 
V, Tc is unchanged and always stay near 3.0ns, that indicates when r — Tcu and Tc2i is 
largely different as in region I and V, AS behavior is exhibited and more prominent in the 
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combination of AS and IS. Furthermore, it is seen in Fig. 9(b) that corresponding Qc is near 
the quahty of only AS (The quahty of only AS is shown by the height of the right peak in 
Fig. 5, where the converse coupling is still not added, and is also lined out in Fig. 9(b) by a 
straight line.) It seems that when two characteristic time scales r — Tc12 and rc2i are largely 
different. The combination effect is not so great to enhance the synchronization quality. In 
region III where Tc2i is close to r — rci2, Tc always approximately equals to Tc2i (Tc ~ Tc2i), 
presenting a striking contrast with Region I and V. The fact that Tc increases with Tc2i in 
region III implies IS behavior is prominent and the decisive force in the combination of AS 
and IS. In addition, the quality of the combination synchronization goes over the straight 
line mostly, indicating the combination of IS and AS has produced a better synchronization 
than only AS in most of region III. The maximum efficiency of the combination is obtained 
at about Tc2i = 3ns (r — = 3ns). Regions II and IV are transition regions. Tc falls 
rapidly from r — r^u to Tc2i in region II, and from Tc2i back to r — r^u in region IV. In this 
transition regions the situation that neither AS nor IS is prominent in the combination may 
occur, where the quality of the resulting combination synchronization is even worse than 
only AS, especially in region II. 



IV. DISCUSSION 



We note that In Fig. 7 where the combination synchronization is demonstrated, the both 
couplings are of the same strength, i.e. 7712 = f]2i- If the feedback in laser 1 is removed, 
the system will turn to be a typical Face-to-Face (F2F) model 0, Q- In F2F model, 
synchronization can be obtained with a time lag between two lasers. And the leader role 
switches from one laser to the other randomly and continuously because of the symmetry. 
Without the explanation from the viewpoint of synchronization combination, it is very hard 
to understand why laser2 will always leads laser 1 after a delayed feedback is added in laser 1. 

A simple analysis of the rate equations for the unidirectional lasers system has identified 
the required matching condition 7712 = 7 (no feedback in laser2) necessary for the existence of 
AS. We note this condition is necessary for the complete synchronization. In fact, when the 
stringent matching condition is not satisfied AS still exists, but it is not so good and always 
hidden behind. As shown in Fig. 7, such not so good AS has even become an important 
factor in the synchronization combination. 
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Sivaprakasam et al experimentally demonstrated a very interesting synchronization phe 



nomenon between two bidirectional coupling lasers |41l. |42|. Their setup configuration is the 
same as Fig. 6 of this paper except for Tc12 = Tc2i in their setup. They found experimentally 
that "slave" can predict the future state of the "master" and the corresponding "anticipating 
time" always equals to the coupling retardation time. We thind the experimentally obtained 
synchronization may arise from the combination of synchronization. And in region III of 



Fig. 9 of this paper, = Tc2i is in good agreement with their experimental result. 



[1] L. M. Pecora and T. L. Carroll, Phys. Rev. Lett. 64, 821 (1990). 

[2] H. G. Winful and L. Rahman, Phys. Rev. Lett. 65, 575 (1990). 

[3] R. Roy and K. S. Thornburg, Phys. Rev. Lett. 72, 2009 (1994). 

[4] T. Sugawara, M. Tachikawa, T. Tsukamoto, and T. Shimizu, Phys. Rev. Lett. 72, 3502 (1994). 

[5] G. D. VanWiggeren and R. Roy, Science 279, 1198 (1998). 

[6] J.-P. Goedgebuer, L. Larger, and H. Porte, Phys. Rev. Lett. 80, 2249 (1998). 

[7] S. Sivaprakasam and K. A. Shore, Opt. Lett. 24, 466 (1999). 

[8] L Fischer, Y. Liu, and P. Davis, Phys. Rev. A 62, 011801(R) (2000). 

[9] H. Pujino and J. Ohtsubo, Opt. Lett. 25, 625 (2000). 

[10] V. Ahlers, U. Parlitz, and W. Lauterborn, Phys. Rev. E 58, 7208 (1998). 

[11] C. R. Mirasso, P. Colet, and P. Garcia-Fernandez, IEEE Photon. Electron. Lett. 8, 299 (1996). 

[12] V. Annovazzi-Lodi, S. Donati, and A. Scire, IEEE J. Quantum Electron. 32, 953 (1996). 

[13] T. Heil, I. Fischer, and W. Elsafier, Phys. Rev. A 58, R2672 (1998). 

[14] C. Risch and C. Voumard, J. Appl. Phys. 48, 2083 (1977). 

[15] C. H. Henry and R. F. Kazarinov, IEEE J. Quantum Electron. 22, 294 (1986). 

[16] T. Sano, Phys. Rev. A 50, 2719 (1994). 

[17] I. Fischer, G. H. M. van Tartwijk, A. M. Levine, W. Elsasser, E. Gobel, and D. Lenstra, Phys. 

Rev. Lett. 76, 220 (1996). 

[18] A. Hohl, H. J. C. van der Linden, and R. Roy, Opt. Lett. 20, 2396 (1995). 

[19] D. W. Sukow, J. R. Gardner, and D. Gauthier, Phys. Rev. A 56, R3370 (1997). 

[20] F. Marino, M. Giudici, S. Barland, and S. Balle, Phys. Rev. Lett. 88, 040601 (2002). 

[21] J. M. Buldii, J. Garci'a-Ojalvo, and M. C. Torrent, Phys. Rev. E 69, 046207 (2004). 



11 



[22] I. Wallace, D. Yu, W. Lu, and R. G. Harrison, Phys. Rev. A 63, 013809 (2000). 
[23] H. U. Voss, Phys. Rev. E 61, 5115 (2000). 
[24] C. Masoller, Phys. Rev. Lett. 86, 2782 (2001). 

[25] Y. Liu, Y. Takiguchi, P. Davis, T. Aida, S. Saito, and J. M. Liu, Appl. Phys. Lett. 80, 4306 
(2001). 

[26] S. Tang and J. M. Liu, Phys. Rev. Lett. 90, 194101 (2003). 

[27] A. Locquet, F. Register, M. Sciamanna, P. Megret, and M. Blondel, Phys. Rev. E 64, 045203 
(2001). 

[28] A. Locquet, C. Masoller, P. Mcgrct, and M. Blondel, Opt. Lett. 27, 31 (2002). 
[29] A. Locquet, C. Masoller, and C. R. Mirasso, Phys. Rev. E 65, 056205 (2002). 
[30] E. M. Shahverdiev and K. A. Shore, Phys. Lett. A 295, 217 (2002). 
[31] A. Murakami and J. Ohtsubo, Phys. Rev. A 65, 033826 (2002). 
[32] L V. Koryukin and P. Mandel, Phys. Rev. E 65, 026201 (2002). 
[33] L. Wu and S. Zhu, Phys. Lett. A 315, 101 (2003). 

[34] S. Sivaprakasam, P. S. Spencer, P. Rees, and K. A. Shore, Opt. Lett. 27, 1250 (2002). 
[35] M. Barahona and L. M. Pecora, Phys. Rev. Lett. 89, 054101 (2002). 
[36] F. M. Atay and J. Jost, Phys. Rev. Lett. 92, 144101 (2004). 
[37] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 (2002). 
[38] R. Lang and K. Kobayashi, IEEE J. Quantum Electron. 16, 347 (1980). 
[39] T. Heil, I. Fischer, W. Elsasser, J. Mulct, and C. R. Mirasso, Phys. Rev. Lett. 86, 795 (2001). 
[40] J. Mulct, C. Masoller, and C. R. Mirasso, Phys. Rev. A 65, 063815 (2002). 
[41] S. Sivaprakasam, E. M. Shahverdiev, P. S. Spencer, and K. A. Shore, Phys. Rev. Lett. 87, 
154101 (2001). 

[42] P. Rees, P. S. Spencer, I. Pierce, S. Sivaprakasam, and K. A. Shore, Phys. Rev. A 68, 033818 
(2003). 



12 



